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Abstract 



We solve the continuous one-dimensional Schrodinger equation for the case 
of an inverted nonlinear delta-function potential located at the origin, ob- 
taining the bound state in closed form as a function of the nonlinear expo- 
nent. The bound state probability profile decays exponentially away from 
the origin, with a profile width that increases monotonically with the non- 
linear exponent, becoming an almost completely extended state when this 
approaches two. At an exponent value of two, the bound state suffers a 
discontinuous change to a delta-like profile. Further increase of the expo- 
nent increases again the width of the probability profile, although the bound 
state is proven to be stable only for exponents below two. The transmission 
of plane waves across the nonlinear delta potential increases monotonically 
with the nonlinearity exponent and is insensitive to the sign of its opacity. 
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The delta-function potential 5{x — xq) has become a familiar sight in the 
landscape of most elementary courses on quantum mechanics, where it serves 
to illustrate the basic techniques in simple form. As a physical model, it has 
been used to represent a localized potential whose energy scale is greater than 
any other in the problem at hand and whose spatial extension is smaller 
than other relevant length scales of the problem. Arrays of delta-function 
potentials have been used to illustrate Bloch's theorem in sohd state physics 
and also in optics, where in the scalar approximation, wave propagation in a 
periodic medium resembles the dynamics of an electron in a crystal lattice. It 
is well known that the single "inverted" delta-function potential — Q S{x — xo) 
possesses one exponentially localized bound state for all values of the opacity 
parameter H,. Its existence and stability has been tested against the effects 
of different boundary conditions^ and symmetry- breaking perturbations^. In 
addition, the inverted delta potential has been used as a semi-permeable 
barrier to examine resonance phenomena in scattering theory^, among others. 
Other interesting applications of the delta-function potential concept are 
found in Ref.4. 

In this work we examine the problem of finding the bound state and the 
transmission coefficient of plane waves across an inverted nonlineair delta- 
function potential, described by the so-called Nonlinear Schrodinger (NLS) 
equation: 



where > is the opacity coefficient and a is the nonlinearity exponent. 




(1) 
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For a = we recover the familiar problem of the linear "inverted" delta- 
potential which possesses an exponentially decaying bound state profile for 
any opacity strength: = exp(— |a;|). The rate of decay in 

space is determined by the localization length 2/VL which increases (decreases) 
as the opacity decreases (increases). 

It might seem odd at first to see a nonlinear-looking SchTodinger equation like 
(|l]), since we know that quantum mechanics is linear. Therefore, all physical 
systems should be described by coupled sets of linear equations. However 
we oftentimes can only concentrate on a few "relevant" degrees of freedom, 
making suitable approximations to deal with the rest. At times, the price to 
pay for this reduction is the appearance of nonlinear evolution equations for 
the variables of interest, such as Eq. (|^). 

For instance, in atomic physics, a well-known approximation when dealing 
with multi-electron atoms is the self-consistent field approximation (Hartree- 
Fock). In this case each electron is described by a single-particle wave func- 
tion that solves a Schrodinger-like equation. The potential appearing in this 
equation is that generated by the average motion of all the other electrons, 
and so depends on their single-particle wave functions. This results in a set of 
nonlinear eigenvalue equations^. A more recent application of the mean-field 
ideas to a weakly interacting Bose condensate can be found in Ref. [6]. 
For a = 2, Eq.(|l|) could model the problem of an electron propagating in a 
one-dimensional linear medium which contains a vibrational "impurity" at 
the origin that can couple strongly to the electron. In the approximation, 
where one considers the vibrations completely "enslaved" to the electron, one 
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obtains Eq.(|l|) as the effective equation for the electron. 
A closely related equation, given by the discrete version of (|^) is known as 
the discrete nonlinear Schrodinger (DNLS) equation. It was introduced in 
its time-independent form, in the late fifties by Holstein in his studies of 
the polaron problem^ in condensed matter physics. The DNLS equation 
was derived in a fully time-dependent form by Davydov in his studies of 
energy transfer in proteins and other biological materials^. In the continuum 
limit the time-dependent DNLS equation reduces to the time-dependent NLS 
equation, which supports soliton solutions. Therefore, a soliton-based energy 
transport appears as an attractive candidate mechanism for energy transport 
in biomolecules. A recent review of the status of Davydov's proposal can be 
found in Ref. [9]. The time-dependent DNLS equation can also be viewed 
as the evolution equation for a Hamiltonian system of classical anharmonic 
oscillators^^. 

An important application of the continuous model (|T]) is that of a wave 
propagating in a one-dimensional linear medium which contains a narrow 
strip of nonlinear (general Kerr-type) material^ ^. This nonlinear strip is as- 
sumed to be much smaller than the typical wavelength. In fact, periodic and 
quasiperiodic arrays of nonlinear strips have been considered by a straight- 
forward generalization of Eq.(|I]) in order to model wave propagation in some 
nonlinear superlattices^^. 

Bound State (E = E^'^^ < 0): Our system consists of a single, infinitely 
localized potential well in a continuous infinite line and therefore, lacks any 
natural length scale. If the delta potential were confined between two infinite 
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walls, the distance between the walls would provide a length scale. If, instead 
of a continuous line, the potential were defined on a discrete lattice, its lattice 
constant would define a natural length scale. Also, if instead of one delta 
potential, we had at least two of them, their mutual distance would constitute 
a natural length scale for the system. 

In our case we have none of these. Thus, Q serves only to define the unit of 
distance (as it does in the linear case). It is possible to get rid of Q formally 
as follows: From Eq.(|l]) we see that f2 must have units of [distance] 
which suggests the definition of a dimensionless distance u as u = x/L, with 
L = In terms of u and (j){u) = (f){x), Eq.(|l]) can be recast 

in a dimensionless form: 



with = -2mL^E^^^/h'^ and = {(P / du^)(t){u). The opacity Vt has now 

disappeared from view since it only determines the unit of distance. We try: 



Using the continuity of and the discontinuity of (j)'{u) at m = 0, one ob- 
tains A = B and k = (1/2) |v4|". Finally, use of the normalization condition 



(P"{u) - e 0(m) 



6{u) i0(«)r 0(«), 



(2) 




M < 
M > 



(3) 



i = i-oo\m\ 



^ du, leads to 




u 



(4) 



with a dimensionless bound state energy 




(5) 
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As in the linear (a = 0) case, the bound state profile is exponentially de- 
creasing away from the delta potential with localization length 2^/*^^""^ As 
a increases from zero, the probability profile widens and the bound state 
energy decreases in magnitude. At a = 2~, the state is completely extended 
all over the real axis and the bound state energy is vanishingly small. At 
a = 2+, the bound state becomes infinitely localized, with a delta-like prob- 
ability profile and with an infinite bound state energy. Further increase in 
the nonlinear exponent leads to a widening of the probability profile and to 
a corresponding reduction in the magnitude of the bound state energy. 
Figure 1 shows the amplitude, or the inverse square probability profile width, 
of the bound state as a function of the nonlinearity exponent. However, at 
this point, an important observation is in order. The total energy of the 
system does not coincide with the bound state energy. In order to see this, we 
must consider the full time-dependent nonlinear Schrodinger equation that 
gives rise to Eq. (|^). By using r = t/T as a dimensionless time variable, 
with T = [Tijlvnl?'^'^ ^ we have 

I — = -— |^(M,r)| V^(M,r). (6) 

(XT a T'^ 

In other words, we have i{d/dT) tp{u, t) = H ip{u, r), where the Hamiltonian 
operator can be decomposed as H = Hq + O^l, where Hq = p^, with p = 
i (d/du) and Onl = —^u) \iIj{u,t)\°' as the nonlinear part. We see that H 
depends on time explicitly, through the time dependence of O^l: 

dr dr ^ ' 
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This implies that {H) is no longer a constant of the motion: 

For the nonlinear part, we have 

But, 



\ J \ I I \\2 



du\^l,{u,T)\\dONL/dT). (10) 



\ Or I 

By expressing Oat^ in terms of ipiu, r) and using Eq. (^, we can recast (^) 
as 

Id0.,\ ^^H,o.^^^, (11) 



dt 2 



which means 



« = i (l + f)(|i/,0„J). (12) 
By comparing Eqs. (|lT]) and ([12|), we conclude 

Finally, by inserting this back into Eq.(|^, Eq.(P) becomes: 



'//-f^)O„i)=0. (15) 



which implies 

dr \" \a + 2 

Therefore, the true energy operator for our system is Ht = H—{a/a+2) Onl- 
For a stationary-state, il){u^ r) = (j){u) exp{—iE^''H/h), the total dimensionless 



energy is then 

Thus, for a < 2 the total energy is negative and the eigenstate is a stable 
localized state. On the contrary, when a > 2, the total energy is positive and 
the eigenstate is localized but possibly unstable, which means that any weak 
'perturbation' could make it disappear into the continuum. This explains the 
'stable' and 'unstable' labelling in Fig.l. Only for a = 0, i.e., the linear case, 
both the total energy and the energy eigenvalue coincide. Figure 2 shows 
some probability profiles for several different values of the nonlinear exponent 
that give rise to true (stable) bound states. This distinction between the 
eigenenergy and the total energy must always be kept in mind when dealing 
with effectively nonlinear systems. 

Transmission of plane waves {E > 0): We now cast Eq.(|T]) as 

iIj"{x) + k'^ ij{x) = -n 6{x) \^{x)\'' ^{x) (17) 

where k'^ = 2mE/li? is the electron wavevector. Unlike the bound state 
problem, we now have 1/k as a natural length scale and can therefore con- 
sider as a bona fide opacity coefficient. The problem looks similar to the 
usual single delta-barrier problem, with the exception of the nonlinear term 
l^/^l" that modulates the strength of the barrier opacity, depending on how 
much electronic probability is sitting on the barrier. We will examine the 
dependence of the transmission coefficient on VL and a. 
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Since we are interested in plane wave transmission, we set 

,/ ^ I Roexp(ikx) + Rexr)(—ikx) x<0 
ip[x] = < „ / ' ^ (18) 

^ ( Texp{ikx) X > ^ ' 

From the continuity of iIj{x) and discontinuity of iIj'{x) at x = 0, we obtain 

T = Ro + R (19) 
ikT = ik{Ro - R) - n |T|" T. (20) 

From here, one obtains T = 2 Ro/{2 — (iVt/k) |T|°). Defining the trans- 
mission coefficient as t = |Tp/|i?oP, we obtain the following equation for the 
transmission coefficient: 

4 

^ " A + {n*/ky ^^^^ 

where Q* = fi|-Ro|" is the "effective" opacity. We note that ( pT] ) is invari- 
ant under a sign change in Q. In other words, both the "upright" and the 
"inverted" delta potentials possess identical transmissivities. 

For arbitrary a, Eq. (^) is a nonlinear equation for t and must be solved 
numerically. There are, however, four exactly solvable cases, three of which 
can be described shortly: 

1. a = (linear case): From ( PT| ) we immediately obtain the well-known 
result 

^ = i + iQ/2ky ^^^^ 

2. a = 1: Now Eq. (|l|) can be recast as the quadratic equation (fi*/2fc)^ 



t — 1 = 0, with physical solution 



(23) 
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3. a = 2: Now we deal with a cubic equation for t: (f2*/2A;)^t^+t — 1 = 0. 
Its physical solution is 

t = (2/9^/^ A{k, n*) - (32/3)^/3 (A;V|^^*|) A{k, (24) 

where A{k,n*) = 9k'^\n*\ + ^3{16k^ + 27k^n*^) 

The case a = 3 is exactly solvable in principle, but it leads to a cumbersome 
expression for t that is not particularly illuminating. 

If we recast the general equation for t as t (1 + (fi*/2fc)^ t") = 1, a bit 
of simple analysis will convince the reader that the left hand side is always a 
monotonically increasing function of t for a > 0. Therefore, there is always 
only one solution in the interval < t < 1. Figure 3 shows the transmission 
coefficient t as a function of A;/r2* and several different nonlinearity exponents 
a. Unlike the bound state calculation, there is no restriction here on the 
magnitude of the nonlinear exponent a. For all wavevectors, the transmission 
increases with increasing a and does not display any special behavior at 
a = 2. The increase of t with a can be easily understood with the help of 
Eq. (|21|) : For any a > 0, < 1 since t is less than unity. Thus, Q* < Q* 
which means that the total "nonlinear" opacity is always smaller than the 
"linear" one, hence a higher transmission. 

Summary. In this work we have calculated the bound state correspond- 
ing to a single "inverted" nonlinear delta-function potential, with opacity Q 
and nonlinearity exponent a. Following the usual methods of elementary 
quantum mechanics, we arrived at a closed form expression for the bound 
state characterized by an exponentially-decreasing probability profile, with 
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a localization length that decreases with increasing a. The most significant 
feature of this solution is the existence of a critical a value, namely 2, beyond 
which the total energy (not the eigenenergy) of the bound state becomes pos- 
itive, making the state unstable against a collapse into the continuum. The 
transmission of plane waves across the nonlinear delta potential is invariant 
under a sign change in opacity, and increases monotonically with an increas- 
ing nonlinearity exponent. The transmission is always higher than in the 
linear case, for a nonzero exponent. 

Finally, it is important to remark that because of the nonlinear nature of 
Eqs. (|l|) and (|), it is no longer possible to superpose stationary states in 
order to find the time evolution of a given initial state. A stationary state 
solution of Eq.(^) is now only a particular solution whose relation to the 
solution of the time-dependent problem is unclear. Other features that arise 
in similar 'nonlinear' quantum mechanical problems include the fact that 
eigenstates are no longer guaranteed to be orthogonal to each other. Also, 
the number of eigenstates is no longer constant, but depends on nonlinearity. 
Thus, 'nonlinear' quantum mechanics is considerably more challenging than 
the linear one, although the reader should be aware that, as was mentioned 
at the beginning of this Note, nonlinearity in quantum mechanics is the 
consequence of some underlying assumption about the system. 
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Captions List 



Fig.l : Normalized bound state amplitude at the origin as a function of the 
nonlinearity exponent. The wavefunction changes discontinuously at a = 2, 
becoming unstable for a > 2. 

Fig. 2 : Bound state probability profile for the "inverted" nonlinear delta- 
function potential, for several nonlinearity exponents a that give rise to a 
stable bound state. 

Fig.3 : Transmission coefficient of plane waves across the nonlinear delta- 
function potential versus wavevector, for several different nonlinearity expo- 
nent values. The transmission is the same for the "upright" and "inverted" 
delta potentials. 
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